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Abstract: The recent Department of Energy Atmospheric Radiation Measurement (ARM) 
Aerosol Intensive Operations Period (AIOP, May 2003) yielded one of the best measurement 
sets obtained to-date to assess our ability to measure the vertical profile of ambient aerosol 
extinction 0^(7.) in the lower troposphere. During one month, a heavily instrumented aircraft 
with well characterized aerosol sampling ability carrying well proven and new aerosol 
instrumentation, devoted most of the 60 available flight hours to flying vertical profiles over the 
heavily instrumented ARM Southern Great Plains (SGP) Climate Research Facility (CRF). This 
allowed us to compare vertical extinction profiles obtained from 6 different instruments: airborne 
Sun photometer (AATS-14), airborne nephelometer/absorption photometer, airborne cavity ring- 
down system, ground-based Raman lidar and 2 ground-based elastic backscatter lidars. We find 
the in-situ measured a ep (7.) to be lower than the AATS-14 derived values. Bias differences are 
0.002 - 0.004 Km' 1 equivalent to 12-17% in the visible, or 45% in the near-infrared. On the other 
hand, we find that with respect to AATS-14, the lidar 0^(7.) are higher: Bias differences are 
0.004 Km" 1 (13%) and 0.007 Km' 1 (24%) for the two elastic back-scatter lidars (MPLNET and 
MPLARM, X=523 nm) and 0.029 Km' 1 (54%) for the Raman lidar (Jl=355 nm). An unnoticed 
loss of sensitivity of the Raman lidar had occurred leading up to AIOP and we expect better 
agreement from the recendy restored system. 

Looking at the collective results from 6 field campaigns conducted since 1996, airborne in 
situ measurements of 0^(7.) tend to be biased slightiy low (17% at visible wavelengths) when 
compared to airborne Sun photometer o cp (X). On the other hand, < 3 fp (X) values derived from 
lidars tend to have no or positive biases. 

From the bias differences we conclude that the typical systematic error associated with 
measuring the tropospheric vertical profile of the ambient aerosol extinction with current state- 
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of-the art instrumentation is 15-20% at visible wavelengths and potentially larger in the UV and 
near-infrared. 

1 Introduction 

A major uncertainty in predicting future changes to die Earth system in general, and its 
climate in particular, stems from the difficulty of modeling the effects of atmospheric aerosols. 
In fact, recent modeling studies debate to what extent controlling the emission of aerosol (i.e. 
reducing the emission of light-absorbing aerosol) into the Earth’s atmosphere may be a feasible 
way to slow global wanning [Jacobson, 2002; Hansen et al., 2000; Sato et aL, 2003; Penner et 
al., 2003; Penner, 2003]. The current low confidence in the estimates of aerosol induced 
perturbations of the Earth’s radiation balance is caused by the highly non-uniform compositional, 
spatial and temporal distribution of tropospheric aerosols owing to their heterogeneous sources 
and short lifetimes. 

Aerosols affect climate through a variety of pathways. These pathways include direct effects 
on the scattering and absorption of radiation, indirect effects caused by aerosol roles in cloud 
microphysics, and “semi-direct” effects caused by aerosol modification of atmospheric heating, 
temperature profiles, convection, and large-scale horizontal transport [e.g., Ackerman et al . , 
2000; Chameides and Bergin, 2002; Lelieveld et al., 2002; Menon et al., 2002]. Many of these 
pathways can affect precipitation, and thus aerosols are intimately linked to the hydrological 
cycle [e.g., Ramanathan et al., 2001; Rotstayn and Lohmann, 2002]. 

Monitoring the global distribution of aerosols requires the combination of continuous 
observations from satellites, networks of ground-based instruments, and dedicated field 
experiments [ Kaufman et al., 2002]. 
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The globally distributed AErosol RObotic NETwork (AERONET) consisting of -200 Sun- 
and sky-scanning ground-based automated radiometers provides column measurements of 
aerosol optical properties, with up to ten years of observations in some locations [Holben et al., 
2001]. These data are used extensively for the validation of satellite-derived aerosol properties 
[e.g. Diner et al., 2001; Torres et al., 2002; Chu et al., 2003]. In situ measurements of aerosol 
optical properties and composition are made by numerous ground-based networks around the 
world [e.g. Delene and Ogren, 2002; VanCuren, 2003]. Ground-based lidar networks monitoring 
the vertical distribution of aerosols are also emerging [Welton et al., 2001, Ansmann et al., 
2003]. The era of continuous satellite-based observation of the vertical distribution of 
tropospheric aerosols has begun very recently with the launch of the Geoscience Laser Altimeter 
System (GLAS) in January 2003 [Spinhime et al., 2003]. 

Here, we assess the accuracy with which the vertical profile of aerosol extinction (a 
fundamental aerosol property) can currently be measured with state-of-the art instrumentation. 
We cannot stress enough that for climate considerations it is the properties of the unaltered 
aerosol at its ambient concentration and thermodynamic state that are of interest. Hence the 
accuracy assessment presented here applies to the measurement of the vertical profile of ambient 
aerosol extinction. To arrive at this assessment we rely on comparisons of ambient aerosol 
extinction profiles obtained in coordinated field campaigns that include in situ and remote 
sensing measurements of aerosols aboard airborne platforms over surface-based Ildars. We start 
with the results of a recent campaign, the Department of Energy Atmospheric Radiation 
Measurement (ARM) Aerosol Intensive Operations Period (AIOP, May 2003), and then consider 
these results in the context of findings from other field campaigns conducted since 1996. 


AIOP yielded one of the best-suited measurement sets obtained to-date to assess our ability 
to measure the vertical profile of ambient aerosol extinction. During one month, a heavily 
instrumented aircraft with well characterized aerosol sampling ability carrying a combination of 
well proven and new aerosol instrumentation, devoted most of the 60 available flight hours to 
flying vertical profiles over the heavily instrumented ARM Southern Great Plains (SGP) Climate 
Research Facility (CRF) [ Ackerman and Stokes, 2003]. This allows us to compare vertical 
extinction profiles obtained from 6 different instruments: airborne Sun photometer, airborne 
nephelometer/absorption photometer, airborne cavity ring-down system, ground-based Raman 
lidar and 2 ground-based elastic backscatter lidars. 

2 Measurements 

2.1 Airborne Measurements 

2.1.1 The Twin Otter aircraft 

The Twin Otter is operated by the Marina, California, based Center for Interdisciplinary 
Remotely-Piloted Aircraft Studies (CIRPAS) [Bluth et al., 1996; Bane et al., 2004]. Between 
May 6 and May 29, 2003, the Twin Otter performed 16 research flights out of Ponca City, 
Oklahoma, Airport. All flight patterns were anchored at the ARM SGP CRF (36.60°N, 97.48°E, 
319 m), 32 km west of Ponca City. For the AIOP campaign the maximum flight altitude was 5.6 
km. 

2.1.2 Aerosol Extinction from Sun photometry aboard the Twin Otter 

The NASA Ames Airborne Tracking 14-channel Sun photometer (AATS-14) measures the 
transmission of the direct solar beam in 14 spectral channels (354 to 2139 nm). AATS-14 is an 
enhanced version of the AATS-6 instrument [Matsumoto et al., 1987]. 
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The AATS-14 tracking head is mounted outside the aircraft skin to minimize blockage by 
aircraft structures and to avoid data contamination by aircraft-window effects. The instrument 
locates and tracks the Sun without input from an operator and records data in a self-contained 
data system. Using aircraft-provided data on latitude, longitude and ambient static pressure, 
aerosol (or particulate) optical depth x p (X) and columnar water vapor (CWV) are computed and 
displayed in real-time. 

AATS-14 made its first science flights during the Tropospheric Aerosol Radiative Forcing 
Observational Experiment (TARFOX) in July 1996 [ Russell et al., 1999a,b]. Since then, AATS- 
14 has been operated on many aircraft in numerous aerosol oriented field experiments: ACE-2 
[Schmid et al., 2000], SAFARI 2000 [Schmid et al., 2003a], ACE-Asia [Schmid et al., 2003b], 
CLAMS [Redemann et al., 2005], SOLVE-2 [Livingston et al., 2005; Russell et al., 2004a]), and 
ADAM [Bucholtz et al., 2003]. 

During AIOP, AATS-14 operated successfully on all 16 Twin Otter research flights. 
Conditions in the boundary layer tended to be relatively turbulent, resulting in larger (compared 
to flights over the ocean surface) AATS-14 tracking errors. Measurements exceeding a tracking 
error of 1° were flagged as questionable data points and not used for this study. The tracking 
capabilities of AATS-14 under such bumpy conditions have recently been improved by changing 
settings in the tracking software. To avoid contamination of the AATS-14 entrance window, the 
tracking head was moved into its park position before flying through clouds. 

Our methods for data reduction, calibration, and error analysis have been described 
previously [Russell et al., 1993a; Schmid and Wehrli, 1995; Schmid et al., 1998 and 2001]. A 
brief s umm ary is given here. The AATS-14 ch ann els are chosen to allow separation of aerosol, 
water vapor, and ozone transmission. From these slant-path transmissions we retrieve t p (L) in 13 
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narrow wavelength bands and the columnar amounts of water vapor and ozone. In addition to the 
corrections for Rayleigh scattering and O 3 absorption, some channels require corrections for 
NO 2 , H 2 O and O 2 -O 2 absorption. Cross-sections were computed using LBLRTM 6.01 [Clough, 
and Iacono , 1995] with the CKD 2.4.1 continuum model using the Hl'l'kAN 2000 (v 1 1.0) line- 
list [Rothman et al., 2001, 2002 ] (including an update for water vapor from 04/2001, see 
http://www.hitran.com/hitran/updates.html). NO 2 cross-sections not included in LBLRTM 6.01 
were taken from Harder et al. [1997]. NO 2 was assumed constant at 2xl0 ' 15 molecules cm' 2 . 

The AIOP AATS-14 dataset consists of 13 wavelengths (354, 380, 453, 499, 519, 604, 675, 
778, 865, 1019, 1241, and 2139 nm) at which we retrieve x p (X.) and the 941-nm wavelength, 
which we use to determine CWV [Schmid et al., 2001]. 

The columnar O 3 content needed to correct for O 3 absorption was derived from high altitude 
(hence low x p (k)) spectra (discussed below) by using the spectral fitting technique introduced by 
King and Byrne [1976] and validated recently by Livingston et al. [2005]. The so-determined 
columnar O 3 content 0 3 (zj) corresponds to the flight altitude z„ at which the low x p (k) spectra 
were measured. Values at all other flight altitudes were determined by scaling a standard O 3 
profile so it passes through 03 (zO- 

AATS-14 was calibrated at the Mauna Loa Observatory (MLO), Hawaii, 1.5 months before 
and 1 .5 months after the AIOP campaign using the Langley plot technique [Schmid and Wehrli, 
1995]. As a result of band-pass filter degradation, the calibration constants obtained from the 
post-mission calibration were slightly different from those obtained from the pre-mission 
calibration. None of the 14 calibration constants had changed by more than 1.6% with 5 channels 
exhibiting a change of less than 0.5%. 
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To determine the best calibration constants, V 0 (X), applicable to the AIOP data set we 
inspected spectra with low t p (E) values measured during higher altitude legs. This resulted in 16 
spectra taken during 14 flights with x p (7.) between 0.01 and 0.06 (at 499 nm) at altitudes 3. 1-5.6 
km. Starting with calibration constants obtained by linearly interpolating Vo(X) between pre- and 
post-mission calibration, we then adjusted the calibration constants within the bounds of pre- and 
post-mission calibration in such a fashion that the retrieved x p (A.) yielded “smooth” x p (7.) spectra 
for all 16 high-altitude cases. This procedure revealed that it is best to use slightly different 
calibration constants for different periods with the AIOP period. This fine-tuning of the 
calibration constants indicates that some of the optical filters must have degraded in a stepwise 
fashion. 

During AIOP, AATS-14 sampled at 3 Hz with data recorded every 4 seconds consisting of 
an average of 9' samples taken in the first 3 of the 4 seconds. The sample standard deviation of all 
science detector outputs is also stored in the data files. These standard deviations were used in 
our cloud-screening algorithm that is based on clouds exhibiting higher standard deviations than 
clear sky. This cloud-screening method can be ambiguous when thick and highly variable dust 
layers are present above the aircraft. However we did not encounter such conditions during 
AIOP. 

Because Sun photometers have a nonzero field of view (FOV), they measure some diffuse 
light in addition to the direct solar beam. As a result, uncorrected Sun photometer measurements 
can overestimate direct-beam transmission and hence underestimate t p (X). For most aerosol 
conditions and Sun photometer FOVs these effects are negligible. For example, Eck et al. [1999] 
report that for the AERONET sun/sky radiometers, which have FOV half-angle 0.6°, the diffuse- 
light correction to apparent x p (A.) is <0.7% of x p (?i), even for desert dust with aerosol effective 
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radius as large as 1.75 [ini. AATS-6 and -14, are designed and built with a relatively large FOV 
(measured half-angle 1.85°) to help keep the full solar disk in view when sun-tracking during 
aircraft maneuvers. This larger FOV makes it necessary to assess quantitatively the diffuse light 
effects on AATS-derived x p (k) when large particles are dominant. We have previously done this 
for postvolcanic stratospheric aerosols [Russell et al., 1993a,b] and for the mineral dust 
dominated Puerto Rico Dust Experiment (PRIDE) [Livingston et al., 2003] and ACE-Asia 
[Schmid et al. 2003b, Redemann et al., 2003] campaigns. Russell et al. [2004a] established 
correction factors that correlate well with aerosol effective radius and also with Angstrom 
exponent 


Ox(A.i ,2.2>=-ln[Tp(Xi )/x p (A.2)]/ln(Xi A-2>, ( 1 ) 

We find the correction factors to be negligible for the Angstrom exponents encountered 
during AIOP. 

The total uncertainty bx p (X) of the retrieved T p (k), due to uncertainties in calibration, sun- 
tracking, signal measurement, airmass computation, and corrections of molecular scattering and 
absorption, was computed following the procedures given by Russell et al. [1993a] and Schmid 
et al [1997]. Note that the impact of tracking errors can be treated as calibration errors taking 
into account the tracking deviation from the Sun and the measured dependence of each channel’s 
response on this deviation angle. In most instances, 5 t p (X) is dominated by the uncertainty in 
Vo(A). Neglecting for the moment the dependence of dz p (k) on the other factors mentioned 
above, one obtains [Russell et al. ,1993a] 


*,(*) = 


1 8V 0 {X) 

m V 0 {A) 


( 2 ) 
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with 


cos 0 

Hence a relative uncertainty of 1% in the calibration constant Vo will lead to an absolute 
uncertainty in the aerosol optical depth Sz p (X) of 0.01 for a solar zenith angle 0=0° and to 

smaller uncertainties at larger 0 . The 5t p (X) values obtained using all uncertainties mentioned 
are part of the archived AATS-14 AIOP data. For the data subset and the wavelengths used most 
prominently in this study, this resulted in average dz p (A = 453, 519, 675, 1558 nm)= 0.003, 

0.004, 0.005, 0.003. Note that 5z p {X) is independent of r p {X) except for diffuse light errors 

which we neglect for this study [Russell et at. 1993a]. 

The uncertainty in CWV was computed following Schmid et al. [1996]. For the subset used 

here this resulted in average 5CWV=0.1 1 g/cm . 

During AIOP the Twin Otter was able to fly as low as 90 m above the land surface, thus 
allowing measurement of virtually the entire overlying atmospheric column. Flying at different 
altitudes over a fixed location allows derivation of layer t p (A.) and layer water vapor LWV. 
Differentiation of x p (X) or CWV data obtained in vertical profiles allows derivation of spectral 
aerosol extinction 0^(1) and water vapor density p w (see section 3.2). 

Because most of die errors in x p (X) or CWV are of systematic nature, they cancel out when 
differences (such as layer x p (X) or LWV) or differentiations (cs rp (X) or p w ) are used. However 
since the aircraft requires a finite time to fly a vertical profile which has a finite horizontal 
component, temporal and horizontal variation of the aerosol above the aircraft will lead to 
uncertainties in the differentiated quantities. The average horizontal variability during AIOP was 
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investigated from t p (X) or CWV measurements during 14 low level legs. This average 
variability, together with the horizontal displacement found during the vertical profiles, was used 
to estimate the uncertainties in the differentiated quantities following the formulas in Rederrumn 
et cd. [2003]. This resulted in average 5a ^(A = 453, 519, 675, 1558 nm)= 0.032, 0.029, 0.024, 

0.014 Km' 1 . 

2.1.3 Aerosol Extinction from Scattering and Absorption Measurements aboard the Twin 
Otter 

Light-scattering data were obtained from four integrating nephelometers aboard the Twin 
Otter. One of these was a three wavelength (450, 550, 700 nm) integrating nephelometer (model 
3563, TSI St. Paul, MN). The other three were Radiance Research (RR) single wavelength (540 
nm) nephelometers (Model RR903, Radiance Research, Seattle, WA). All four were calibrated 
against particle-free air and CO 2 before and at multiple times during the field deployment and 
were zeroed with particle-free air before each flight. All of the nephelometers sampled from a 
shrouded intake whose nominal 50% cutoff diameter was determined to be 8 pm (determined by 
comparison of cross-calibrated interior and exterior FSSP-100 optical probes [Gao et al., 2003]). 

The TSI nephelometer was operated at a flow rate of 30 1/m and with its inlet heater 
operational at -35 °C. This resulted in the RH inside the instrument being considerably lower 
than the ambient RH. The RH inside the TSI nephelometer ranged from near 0 to 35% 
depending on ambient RH. 

The hygroscopic behavior of the aerosol was determined from the three RR nephelometers 
operating at different RH. The three RR nephelometers were operated at RHs below ambient, 
near 85% and at an intermediate level at a flow rate of 6 1/m. The dependence of light-scattering 
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on RH, was parameterized by the exponent of equation (4), based on the work of Kasten [1969] 
(see also Gasso et al., [2000]). 


a sp (RH) = a sp (RH 0 ) 


f 100 -RH v 
100 -RH 0J 


(4) 


where the zero subscript refers to some low, reference RH, and the exponent, y, for the 
measured dependence of light-scattering on RH, is determined by fitting the data to equation (4) 
as in Gasso et al. [2000]. 

We then utilized y to correct the low RH TSI nephelometer scattering signals to the 
measured ambient RH. Though strictly, the determined y would apply only to the wavelength of 
the RR nephelometers (540 nm), we applied it to all 3 TSI nephelometer wavelengths. 

Prior to the humidification correction, the TSI nephelometer a sp (A) values were corrected for 
angular truncation and non-lambertian illumination based on the Angstrom exponent 

aspau^-lnfcsptf.i)/ a sp (2. 2 )]/ln(A.iA.2), (5) 

as suggested by Anderson and Ogren [1998]. No equivalent correction was done for the RR 
nephelometers because their truncation parameters have not been determined. However, since 
the humidograph data are used in a relative sense and dominated by submicrometric particles this 
is not a large error. 

Aerosol light absorption crap(A.) was measured using an improved version of the 3- 
wavelength filter-based absorption photometer (A. = 467, 530, 660 nm) described by Virkkula et 
al. [2005]. The data reduction and correction scheme of Bond et al. [1999] was applied. Because 
CTap(A) was measured just downstream of the TSI nephelometer, it was measured under sub- 
ambient RH (i.e. the same RH as inside the TSI nephelometer to minimize RH dependent 
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artifacts due to the filter substrate). However, following Hegg et al. [1997], no correction was 
made for the higher RH of the ambient air since experimental data for such a correction are 
lacking. A study modeling sulfates with black carbon cores by Redemann et al. [2001] suggests 
that absorption humidification factors are negligible for a wide range of atmospheric conditions. 
However, this may not apply to the considerably more complex real-world aerosol. 

The resulting 0^(7.) and the nephelometer cr^X) were adjusted from temperature Tj and 
pressure p* inside the instruments to ambient (outside the aircraft) T amh and pan* by multiplying 
them with the factor 

k= Pomb_.J±_ ( 6 ) 

Pi '^'amb 

Because the cabin of the Twin Otter is not pressurized, pi is only slightly higher than pan*, 
however Tj is always larger than Tan*. 

The reported nephelometer a sp (X) values were adjusted from their blue, green and red center 
wavelengths (450, 550, 700 ran) to those of the PSAP instrument (467, 530, 660 ran) using the 
Angstrom relationship in Eq. (5). For the comparisons shown in this study, the PSAP crap(^) and 
nephelometer <Tsp(A) were adjusted separately (again using an Angstrom relationship) to 453, 519 
and 675 nm to match AATS-14 and Cadenza (see next section) wavelengths. Aerosol extinction 
was then calculated as 

Oep( A.) = <Tsp(A.) + CTap(k) C 7 ) 

2.1.4 Aerosol Extinction from Cavity-Ring-Down Measurements aboard the Twin Otter 

First demonstrated by O’Keefe and Deacon [1988], the cavity ring-down (CRD) technique 
has been used pr imari ly for gaseous absorption spectroscopy (see various papers in Busch and 


13 



Busch [1999]). The use of CRD to measure aerosol extinction is relatively new [ Smith and 
Atkinson, 2001, Strawa et al. 2003]. The principle behind CRD can be best described using the 
so-called ‘ping-pong’ model. A pulse of laser light is injected into a cavity that consists of two 
highly reflective mirrors. The mirror reflectivity is typically better than 99.96%. The laser pulse 
bounces between the two mirrors inside the ring-down cavity like a ping-pong ball. Each time 
the pulse interacts with the back mirror, a small amount of light (e.g., 0.04%) leaks out. This 
light is collected and detected with a photomultiplier or similar detector. The intensity of the 
light leaking out of the back of the ring-down cavity decreases exponentially. It can be shown 
that the exponential decay, or ring-down time, is related to the mirror reflectivity and the 
extinction of the material inside the cavity. The extinction coefficient is then obtained by the 
difference between measurements made when the cell contains filtered air and when the cell 
contains a particulate-laden flow: 




( 8 ) 


where c is the speed of light, and tc p and kq are the ring-down times of the aerosol laden flow 


and filtered air, respectively. 

Cadenza is the first airborne CRD instrument able to measure aerosol optical properties. The 
prototype Cadenza instrument as described by Strawa et al. [2003] participated successfully in 
the Reno Aerosol Optics Study (RAOS) [Sheridan et al., 2005]. Cadenza then flew its first and 
second successful airborne missions in the ADAM and AIOP experiments aboard the CIRPAS 
Twin Otter. Detailed descriptions of the instrument, the data analysis and comparisons with other 
methods during AIOP are reported by Strawa et al. [this issue]. 
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Using the same aerosol inlet as the nephelometers, Cadenza operated successfully on all 16 
AIOP science flights continuously measuring at X=675 and 1550 nm. Cadenza also measures 
the aerosol scattering coefficient at A.=675 nm. The scattering measurements are discussed by 
Strawa et al. [this issue]. For one minute out of every six minutes Cadenza sampled filtered air. 
The so-derived Ko values were then linearly interpolated to the times when particle-laden air was 
sampled and o ep (X) is determined according to Eq. (8). While not deliberately heated, the sample 
air inside Cadenza was nearly at the temperature of the aircraft cabin and consequently drier than 
the ambient air. Part of this was caused by ram heating at the aerosol inlet and part was due to 
heating of the sample line as it carried aerosol from the inlet to the instrument. We then utilized y 
along with Eq. (4) to correct the low-RH Cadenza u ep (k) to outside-the-aircraft RH and also 
applied the factor in Eq. (6) to correct to outside-the-aircraft pressure and temperature. Though, 
strictly, the y was determined from scattering measurements with the RR nephelometers (7^=550 
nm), we applied it to the Cadenza measurements at h = 675 nm (a^ and scattering portion of a ep ) 
and to 0^(1550 nm). 

2.1.5 Routine Small Aircraft in situ Measurements 

Since March 2000, ARM has been measuring in situ aerosol profiles (IAP) by performing 
routine flights (2-3 times per week) with a small aircraft (Cessna C-172N) over the SGP site. The 
aerosol instr ume nt package consists of a 3-wavelength TSI nephelometer and a PSAP both 
measuring at low RH. There is a 1-pm impactor upstream of the aerosol instruments 
corresponding to a geometric size cut of approximately 0.79 pm (for a particle density of 1.6 
g/cm 3 ). Although the IAP project was not designed to measure ambient a ep (X), Andrews et al. 
[2004] have applied (altitude-independent) corrections for low RH, impactor loss, and limited 
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aircraft ceiling (using information from ground-based nephelometers and Raman lidar) to 
compare the column-integrated LAP <r e? (A.) to the x p (k) measured by ground-based Sun 
photometers. They find the LAP x p (550 nm) to have a consistent offset of -0.04. 

During AIOP the Cessna flew 14 of its standard flights (i.e. level legs at 9 altitudes between 
467 and 3660 m). During five of these flights, the Twin Otter trailed the Cessna on its standard 
legs. This allowed for detailed inter-aircraft comparisons which are presented in companion 
papers by Hallar et al. [this issue] and Andrews et al. [this issue]. 

2.2 Ground-Based Measurements 

2.2.1 Sun photometers 

Three ground-based Sun photometers were used to validate AATS-14 x p (A.) during low 
altitude flybys and to constrain elastic backscatter lidar retrievals. Two of the Sun photometers 
were AERONET Sun and sky-scanning instruments [Holben et al., 1998, 2001; Eck et al., 1999]. 
One of the AERONET instruments (#98) is a standard Cimel CE-318 instrument (providing 
x p (X) at X=340, 380, 440, 500, 670, 870 and 1020 nm) that is operated continuously at the SGP 
CRF. Its data are cloud-screened and quality controlled according to Smirnov et al. [2000]. The 
other AERONET instr ume nt (#125), an extended-wavelength prototype version with an 
additional channel at A^=1640 nm, was deployed specifically for the AIOP. An updated 
processing scheme was applied to the data from AERONET instrument #125 [ Smirnov , 2004]. 

The third Sun photometer was a Normal Incidence Multi-Filter Radiometer (NIMFR). The 
instrument consists of a Multi-Filter Rotating Shadowband Radiometer (MFRSR, Harrison et 
al., [1994]) “head” to which a collimating tube (FOV 5°) is attached. The NIMFR is mounted on 
a solar tracker. z p (k) at five wavelengths (k = 415, 500, 615, 673 and 870 nm) are reported every 
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20 sec. The data set is cloud screened rigorously based on the stability of x p (A.) over about a 10- 
min period using stability limits that were scaled according to the magnitude of x p (A.). 

2.2.2 Micropulse Lidars 

The Micro-Pulse Lidar (MPL) [Spinhime et al., 1995; Campbell et al., 2002] is a single 
channel (X =523 nm), autonomous, eye-safe lidar system originally developed at the NASA 
Goddard Space Right Center and now commercially available. One of the MPLs (hereafter 
referred to as MPLARM) is permanently deployed at the ARM SGP CRF. The second MPL was 
deployed in support of AIOP as part of the NASA Micro-Pulse Lidar Network (MPLNET) 
[Welton et al., 2001], a network of ground-based MPL systems co-located with AERONET 
sun/sky radiometers. 

Vertical profiles of extinction and backscatter were retrieved independently from both co- 
located MPL systems. The retrieval of independent extinction and backscatter profiles from 
single-wavelength elastic backscatter lidar (such as an MPL) faces an inherently ill-posed 
problem, in that it requires the extraction of two unknowns (extinction and 180°-backscatter 
coefficients) from one measurement (the attenuated 180°-backscatter signal) [Ansmann et al., 
1990; Ackermann, 1998]. However, by assuming a constant value of the extinction-to- 
backscatter ratio (S p ) throughout an aerosol layer, and by constraining the integrated extinction 
profile against an independently determined layer t p (X), it is possible to retrieve a unique 
solution for the extinction and backscatter profiles and calculate a layer-averaged value for S p 
[Welton et al., 2000]. This technique yields reasonable results when the atmosphere is well 
mixed, but may produce over- or underestimates of extinction at a given altitude when aerosol 
properties are highly stratified [Welton et al., 2002], 
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The retrievals of MPLARM and MPLNET assume an altitude-independent extinction-to- 
backscatter ratio, S p . For the total column aerosol optical depth, the MPLARM processing uses 
cloud-screened x p (X) retrieved from the NIMFR (discussed above), while MPLNET processing 
uses similarly screened x p (X) [see Smirnov et al., 2000] from the AERONET Sun/sky radiometer 
located at the ARM SGP CRF (also described above). In as much as the MPL systems, their 
calibration and constraining t p (k) were completely independent, the retrievals from MPLARM 
and MPLNET represent independent determinations using fundamentally similar retrieval 
techniques. 

22.3 Raman Lidar 

The CRF Raman lidar (CARL) measures backscattered light at the laser wavelength of 355 
nm as well as the water vapor and nitrogen Raman shifted returns, at 408 and 387 nm, 
respectively. a cp (355nm) profiles are computed from the derivative of the logarithm of the 
Raman nitrogen signal with respect to range [Ansmann et al., 1990]. Unlike with elastic 
backscatter lidars, the Raman technique allows the derivation of Oep profiles without making an 
assumption about the profile of the lidar ratio, S p , and without using the total column x p as a 
constraint [Ansmann etal., 1990; Ferrare etal., 2001]. 

In April 1997, CARL started to operate at the SGP site as a turnkey, automated system for 
unattended, around-the-clock profiling of water vapor and aerosols. To facilitate data processing, 
algorithms were developed to run autonomously delivering water vapor mixing ratio, RH, 
aerosol scattering ratio, aerosol backscatter coefficient, o t , p , and linear depolarization ratio, as 
well as integrated values CWV and x p [Turner et al., 2001, 2002]. The water vapor measurement 
performance of CARL has been characterized extensively (see references in Ferrare et al. [this 
issue]). However, initial comparisons of x p and o ep have revealed discrepancies among the 
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routine CARL, Sun photometer, and the routine small aircraft in situ measurements described 
above [Ferrare et al., 2003]. AIOP was conducted in part to resolve these discrepancies. 
Unfortunately, a gradual loss of the sensitivity of CARL starting about the end of 2001 went 
unnoticed until after AIOP. In an attempt to reduce or remove these adverse impacts, the 
automated algorithms were modified and the AIOP data were reprocessed. Major modifications 
that were made to CARL in 2004 (after AIOP) have dramatically improved the system’s 
sensitivity. This is discussed in more detail by Ferrare et al. [this issue]. 

3 Results / 

In what follows, we will use the AATS-14 measurement of x p (k) and Oep(k) as a reference 
against which we will compare all other methods. This choice is driven by the fact that AATS-14 
has the largest spectral coverage and can match most of the other instruments’ wavelengths 
relatively closely. 

3.1 Comparing x p (k) Obtained from AATS-14 and Ground-based Sun photometers 

As done in previous airborne campaigns, we assess the in-flight performance of AATS-14 
by comparing against surface based Sun photometers. During most of the flights the Twin Otter 
flew at least one low-altitude leg (-90 m above ground) near the SGP CRF. We compared the 
AATS-14 x p (k) with those from the AERONET and NIMFR instruments. During 18 such low- 
altitude fly-bys the AATS-14 data indicate that the direct beam was not obstructed by clouds. In 
two cases involving instrument #125, and three cases involving instrument #98, the 
corresponding AERONET observations had been screened out in the level 2.0 data. For these 
cases we reverted to the non-cloud screened level 1.0 data. It appears that the NIMFR cloud 
screening is even more conservative, in that only 12 fly-bys had concurrent NIMFR data. The 
results of the x p (k) comparison are shown in Table 1. 
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The level of agreement between AATS-14 and the AERONET instruments is similar to what 
we found from low altitude fly-bys over AERONET sites in previous campaigns (i.e. PRIDE 
[Livingston et al, 2003], SAFARI 2000 [Schmid et al., 2003a], and CLAMS [Redemann et at, 
2005]. The agreement between AATS-14 and NIMFR found in AIOP is particularly good, in fact 
operating four Sun photometers (including AATS-6) side-by-side on the ground in previous 
ARM IOPs did not result in a higher level of agreement [Schmid et al., 1999]. 

3.2 AATS-14 Vertical Profiles 

During AIOP, AATS-14 measured numerous vertical profiles of x p (k) and CWV. After 
discarding profiles influenced by considerable spatial inhomogeneity or overlying clouds, we 
derived spectral aerosol extinction o v (X) for 26 profiles by differentiating the x p (k) profiles. 
CWV can be determined despite thin overlying clouds, resulting in 35 CWV and water vapor 
density (p w ) profiles in AIOP. With very few exceptions, the profiles were located directly above 
the SGP CRF. Figure 1 shows 25 x p (L) vertical profiles. Figure 2 shows the corresponding 
profiles. The profiles of CWV for the same 25 cases and the corresponding p w profiles are 
depicted in Figure 3 and Figure 4. To facilitate comparisons, we plotted all profiles on the same 
scale. Gaps in the x p (X) or CWV vertical profiles are caused by temporary blockage of the direct 
solar beam by aircraft structures (tail, antennas) or clouds. 

Most vertical profiles were acquired within 20 minutes of flight time. Occasionally, x p (k) or 
CWV decreased (increased) when the plane descended (ascended). In a horizontally 
homogeneous, time-invariant atmosphere, this would be impossible. However, in the real 
atmosphere it can occur because (1) the Sun photometer can only measure the transmittance of 
the Sun photometer-to-sun path, (2) that path in general passes through a horizontally 
inhomogeneous, time-varying atmosphere, and (3) the path and the atmosphere move with 
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respect to each other as the aircraft moves and the wind blows. Before the Sun photometer i p (k) 
or the CWV profile is vertically differentiated to obtain a ep (A.) or p w , it has to be smoothed (in a 
non-biased manner) to eliminate increases in i p (k) or CJWV with height. In this study we first 
averaged the t p (X) or CWV values over 20-m altitude bins and then used smoothed spline fits for 
this purpose. However, to avoid over-smoothing at altitudes that exhibit actual variations of t p (A.) 
or CWV we occasionally allow u ep (k) or p w to become slightly negative. This can be seen, for 
example, in Figure 2 (top row, 2 nd panel) and Figure 4 (bottom row, 4 th panel). 

Some of the profiles in Figure 2 show elevated aerosol layers with G ep (k) values exceeding 
those in the boundary layer. On May 9, 2003, the aerosol in the elevated layers originated from 
fires in Mexico [ Wang et al., 2004]. The elevated layers observed from May 25 - May 28, 2003 
can be traced back to Siberian fires [Colarco et al., this issue]. The smoke from the intense 2003 
Siberian biomass burning season ultimately traveled around the globe [Damoah et al., 2004]. 

3.3 Comparison of Water Vapor Profiles 

An aircraft in situ-measurement of p w is more straightforward than measuring ambient 
<y ep (k). Several redundant sensors aboard the Twin Otter measured static temperature T, static 
pressure p, and dewpoint temperature T d . from which we computed p w using an expression given 
by Bogel [1977]. 

Since the same vertical differentiation procedure is used to derive <J ep (k) and p w from the 
columnar data r p (k) and CWV, comparing p w obtained from AATS-14 and the aircraft in situ 
sensors should allow conclusions on the robustness of the AATS-14 differentiated profiles of 
p», and rJep(k). 
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In Figure 4, we compare 25 (of 35) vertical profiles of p w derived from AATS-14 and an 
EdgeTech 137-C3 chilled mirror sensor. We observe excellent correspondence between the two 
measurements. This also demonstrates that the differentiated column method can successfully 
reproduce thin (-500 m) dry or humid layers. Figure 5 shows a scatter plot containing all data 
pairs from all 35 profiles. Figure 6 shows a comparison of the layer water vapor (LWV) 
amounts. LWV is obtained by integrating the in situ measured p w over the vertical span of the 
profile and for AATS-14 by subtracting the CWV measured at the top of the profile from CWV 
measured at the bottom. The complete statistics of the comparison are shown in Table 2. The 
agreement in this study is better than what we found during ACE-Asia [Schmid et al., 2003b] 
using the same instrumentation (i.e. in ACE-Asia the rms difference in p w and LWV was 25% 
and 17% vs. 20% and 7% in AIOP). We attribute this to the fact that the AATS-14 AIOP data 
were acquired using a different brand 941-nm filter which was delivered with potentially more 
accurate spectral band-pass information. 

This study finds the Twin Otter chilled mirror p w to be biased slightly high (5%) with 
respect to AATS-14. More extensive AIOP water vapor comparisons are discussed in the 
companion paper by Ferrare et al. [this issue]. 

3.4 Comparison of Aerosol Extinction Profiles 

For the extinction comparison, the profiles from the six methods were binned in 20-m 
altitude bins between 0 and 8 km above sea level. Naturally, empty bins were excluded from the 
comparisons. In virtually all of the comparisons the AATS-14 values were used as the 
independent variable x, however the linear regressions were established using the linear least 
squares bi-sector (lsq-bs) method which minimizes the quadratic distances to the regression line 
in x and y directions [Sprent and Dolby , 1980]. 
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The Nephelometer+PSAP extinctions were adjusted to the closest AATS-14 wavelengths. In 
contrast, the CARL, MPL and Cadenza instruments’ wavelengths are matched closely enough by 
an AATS-14 wavelength that no further adjustment is required. This resulted in eight 
comparisons between AATS-14 and the other five methods at five different wavelengths (see 
Table 3). 

Plotting the profiles allows a visual evaluation on a profile-by-profile basis. Figure 7 makes 
such a comparison for aep(675nm) from AATS-14 and Nephelometer+PSAP. In this 
representation, the Cadenza a e?J (675nm) profiles are virtually indistinguishable from 
Nephelometer+PSAP data points and are therefore not plotted in Figure 7. The high correlation 
(r 2 =0.%3) between the two in situ measurements is evident in the scatter plot representation in 
Figure 8. Averaged over all profiles, the Cadenza a ep (675nm) are higher by 4.7% (based on lsq- 
bs regression line, slope) or 6.6% (based on bias) than the Nephelometer+PSAP values. Strawa et 
al. [this issue] extend the comparison to all 8-sec averages measured during AIOP (not only the 
26 altitude-binned profiles studied here) and find the Cadenza a ep (675 nm) to be higher by only 
0.8% (based on slope of standard regression line forced through origin). This result is obviously 
a very successful demonstration of the airborne application of the CRD method to measure 

Oepfr)- 

Figure 7 shows cases where AATS-14 a ep (615 nm) are in good agreement with the in situ 
measurements, cases where the AATS-14 values oscillate around the in situ data and cases where 
the AATS-14 values are higher. As an illustration, the scatter plot in Figure 9 shows that 
Cadenza cr ep (675 nm) are 10% (based on lsq-bs slope) to 13% (bias) lower than the AATS-14 
values. 
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An alternative way of assessing potential biases in extinction profiles lies in comparing layer 
Layer x p (X) is obtained by integrating the in situ or lidar measured ct ep(X) over the vertical 
span of the profile and for AATS-14 by subtracting the x p (X) measured at the top from the x p (A.) 
measured at the bottom of the profile. As an example, the scatter plot in Figure 10 shows that 
Cadenza layer x p (675 nm) are lower by 15% (based on lsq-bs slope) to 16% (bias) than the 
AATS-14 values. The layer x p (X) comparisons from all methods are summarized in Table 4. 

The comparisons in Table 3 and Table 4 show that the in situ methods yield consistendy 
lower a ep (X) and layer x p (A.) than AATS-14. All regression lines exhibit slopes smaller than 1 
with very small intercepts indicating a proportional difference rather than a systematic offset. 
Based on the slopes, we find the Nephelometer+PSAP CTq,(A.) to be lower by 7%, 10% and 14% 
(A, = 453, 519, 675 nm). The Cadenza ry ep (k) are lower by 10% (X=675 nm) and 39% (A^=1550 
nm). These slopes, the slopes in the layer AOD comparison and also the relative biases in 0^(1) 
and layer x p (X) show a distinct wavelength dependence: The low bias of the in situ measurement 
with respect to AATS-14 increases with increasing wavelength. Partial loss of larger particles 
during sampling would cause the observed spectral behavior. However, so far we have not 
considered a potential wavelength dependence of the humidification correction in Eq. (4). 
Indeed, a 1-year analysis (March 2000 - February 2001) of surface-based dry and hu mi di fi ed 
Osp(X = 450, 550, 700 nm ) (submicron particles only) measured with a TSI nephelometer at the 
SGP CRF [Sheridan et al., 2001], shows a distinct wavelength dependence [Sivaraman et al., 
2004]. As illustrated in Figure 11, this suggests that the humidification correction we applied to 
<Jsp(X) might be an overcoirection in the blue and an under correction in the green (and 
potentially in the near infrared). Kotchenruther et al. [1999] determined a similar wavelength 
dependence of the humidification factor from airborne measurements during TARFOX. Hence, 
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part of the spectral behavior of the low bias observed in this study could stem from a not entirely 
adequate humidification correction. 

Lidar data concurrent with a Twin Otter vertical profile were available in 1 1 (CARL), 13 
(MPLNET) and 19 (MPLARM) cases. Comparison of CARL data with AATS-14 and in situ 
data is discussed in detail in a companion paper [Ferrare et al., this issue]. As summarized in 
Table 3 and Table 4, CARL 0^(355 nm) and layer r p (355 nm) are significantly higher than the 
AATS-14 values. The lsq-bs regression line between AATS-14 and CARL G ep (355 nm) reveals 
an intercept of 0.024 Km' 1 indicating a systematic offset. The mean difference between the two 
data sets is 0.029 Km ' 1 or 54% for the average cr ep (354 nm) of 0.053 Km' 1 . We believe that this 
high bias was pr imari ly due to the unnoticed loss of sensitivity of CARL leading up to AIOP; 
this reduction in sensitivity led to increased calibration errors, larger random errors, and greater 
uncertainties in maintainin g proper alignment, all of which contributed to these differences. We 
expect better agreement in future comparisons from the recently upgraded CARL system. 

Figure 12 shows vertical profiles of from Nephelometer+PSAP and the two Micro Pulse 
Lidars. The three data sets show good agreement for the vertical distribution of aerosol layers 
including fairly thin layers. However, the absolute magnitudes of cr ep (519/523 nm) differ. As 
shown by Table 3, the lsq-bs regression lines of o ep (523 nm) of MPLNET vs. AATS-14 and 
MPLARM vs. AATS-14 reveal intercepts of 0.005 and 0.011 Km' 1 revealing systematic high 
biases. The bias difference is 0.004 Km' 1 (13%) between AATS-14 and MPLNET and 0.007 
Km' 1 (24%) between AATS-14 and MPLARM for the average 0^(519 nm) of 0.030 Km' 1 . 
Surprisingly, the layer r p (5 19/523 nm) comparisons (Table 4) show high biases with respect to 
AATS-14 of 0.023 (MPLNET) and 0.025 (MPLARM) that exceed the biases between AATS-14 
and AERONET#98 (0.008) and NIMFR (0.006) to which the MPL retrievals are anchored. 
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Figure 13 shows the cumulative integrals (top-to-bottom) of the MPL a ep (523 nm) which is 
equivalent to a x p profile as measured by AATS-14. We find that the cumulative integral of 
MPLNET and MPLARM a ep (523 nm) at the top of the AATS-14 profile average 0.014 lower 
than the AATS-14 x p (519 nm). This discrepancy is likely caused by the fact that the MPL 
retrievals determine a ma ximum layer height Zmax (typically below 8-10 km) above which 
<7^(523 nm) is set to 0. The retrievals then assume that the integrated extinction between the 
surface and Zmax make up the total column x p (523 nm) to which they are anchored. Neglecting the 
cumulative a ep (523 nm) from top-of-the-atmosphere to Zmax will therefore introduce a slight high 
bias in the MPLNET a«p(523 nm) and layer x p (523 nm) retrievals. Figure 13 shows that, apart 
from the small bias discussed above, the x p profiles start out at similar values at the top of the 
AATS-14 profiles and end a similar values at the bottom due to the MPL retrievals’ anchoring. 
In between there are discrepancies, though, indicating that the a ep are distributed differendy over 
the vertical profile. This is apparent in Figure 12 where, compared to Nephelometer +PSAP 
<r ep (523 nm), MPLARM indicates lower a ep (523 nm) in elevated layers above 3 km, but higher 
values below 2 km. This may be an effect of the MPL retrievals assuming that S p is altitude 
independent and/or due to inadequate corrections for overlap or afterpulse. 

4 Results from Previous Campaigns 

4.1 Aerosol Extinction from Scattering and Absorption Measurements 

In n umer ous field cam p ai gns since 1996 we have compared a ep (k) and layer x p (>.) obtained 
from Nephelometer+PSAP measurements with either AATS-6 or -14. As shown in Table 5, 
seven data sets from six field campaigns were reported in eight studies. The data sets were 
obtained aboard five different airplanes. Different metrics have been used in the eight studies to 
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describe the level of agreement making direct comparisons difficult. For each study we have re- 
computed the coefficient of determination, r 2 , and the relative bias differences in layer x p (A.). The 
studies cover wavelengths between 450 and 700 nm. We find that with respect to AATS-14 (or - 
6) the layer x p (A.) value from Nephelometer+PSAP are biased low by -5 to -33%. The average 
low bias (all A.) is -17%. 

Several studies have compared column integrated a ep (A.) with ground-based Sun photometer 
measurements of x p (A.) ( Remer et al. [1997]; Kato et al. [2000]; Andrews et al. 
[2004]). Invariably they find the in-situ derived x p (A.) to be biased low with respect to the Sun 
photometer measurements. However, in all three studies assumptions about the aerosol above the 
maximum aircraft samp ling altitude had to be made, humidification factors were riot measured 
on the aircraft, and aircraft inlets were not suitable for sampling of larger aerosol particles. 

4.2 Extinction calculated from airborne particle size distributions 

In three campaigns since 1996 we have compared a ep (k) and layer x p (A.) calculated from 
airborne measurements of particle size distributions with either AATS-6 or -14. As shown in 
Table 6, excellent agreement was achieved in ACE-2 and ACE-Asia but poorer agreement 
resulted from the TARFOX data set. While there is also a tendency for the o ep (k) and layer x p (A.) 
calculated from in situ data to be lower than the AATS-14 values, we observe the low bias found 
in ACE-2 and ACE-Asia to be smaller than in the corresponding Nephelometer+PSAP 
comparisons listed in Table 5. However, the reverse is the case for the TARFOX da ta set. 

To our knowledge, Clarke et al. [1996] present the only other study where layer 
t p (X) calculated from particle size distributions were compared to the values obtained with an 
airborne sunphotometer (different from AATS). Good agreement was achieved for a profile 
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dominated by pollutant aerosol but the calculated 0^(550 nm) values were 50% lower in a 
profile featuring an elevated Saharan dust layer. 

4.3 Extinction from surface based and airborne lidars 

In numerous campaigns we have compared a £p (A.) and x p (A.) vertical profiles from one of the 
AATS instruments with surface-based or airborne lidars (see Table 7). The results involving 
seven different types of lidar systems have been published in 10 studies. Unfortunately we found 
it difficult to convert the results from all 10 studies into one quantitative metric. Therefore we 
use qualitative terms to describe the bias differences in cs ep (k). As can be seen from Table 7, 
many comparisons result in small or no biases, however the biases that do occur are positive (i.e. 
lidar cr v (X) larger than AATS values). 

Masonis et al. [2002] compared aircraft in situ and Raman lidar profiles of a ep and 180° 
backscattering during the 1999 Indian Ocean Experiment (INDOEX). They found the lidar- 
derived values to be -30% larger than the in situ-derived values. Petzold et al. [2002] report 
agreement wi thin 30% (rms) for 180° backscattering measured with a six wavelength 
Raman/Mie lidar and calculated from airborne size distribution measurements during the 
Lindenberg Aerosol Characterization Experiment (LACE 98). 

5 Summary and Conclusions 

AIOP yielded one of the best measurement sets obtained to-date to assess our ability to 
measure the vertical profile of ambient aerosol extinction a ep (X). Extensive vertical profiling of 
the CIRPAS Twin Otter, carrying state-of-the art aerosol and radiation instrumentation, over the 
heavily instrumented ARM CRF allowed us to compare 1 1 to 26 a ep (h) profiles obtained from 6 
different instruments: airborne Sun photometer (AATS-14), airborne nephelometer plus 
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absorption photometer (Nephelometer+PSAP), airborne cavity ring-down system (Cadenza), 
ground-based Raman lidar (CARL) and 2 ground-based elastic backscatter lidars (MPLARM and 
MPLNET). 

We find good agreement among the in-situ measurements, Cadenza and 
Nephelometer+PSAP, on the Twin Otter aircraft. Averaged over 25 profiles, the Cadenza 
07^(675 nm) are higher by 6.6% (bias difference) than the Nephelometer+PSAP values. This 
represents a very successful demonstration of the first airborne application of the cavity- 
ringdown method to measure <Jep(k). 

Subsequently we used the AATS-14 measurement of t p (A.) and <y ep (X) as a reference against 
which we compared all other methods. This choice was driven by the fact that AATS-14 has the 
largest spectral coverage and can match most of the other instruments’ wavelengths relatively 
closely. 

When compared to AATS-14 a ep (k), we find the in-situ measurements to be biased low 
(0.002 - 0.004 Km' 1 equivalent to 12-17% in the visible, or 44% in the near-infrared). The low 

bias is also apparent when considering layer T p {X). The statistical quantities we investigated 

( 

show that the differences (which should be considered modest, at least for the visible) are 
proportional differences rather than systematic offsets. We also find the low bias to increase with 
increasing wavelength. 

On the other hand, we find that with respect to AATS-14, the o ep (X) values from all 3 lidars 
are biased high: Bias differences are 0.004 Km’ 1 (13%) and 0.007 Km 1 (24%) for the two elastic 
back-scatter lidars (MPLNET and MPLARM, X=523 nm) and 0.029 Km 1 (54%) for the Raman 
lidar (k = 355 nm). Unlike the differences found between AATS-14 and the in-situ measurements, 
the differences between AATS-14 and the three lidars have the nature of an offset. This causes 
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the relative bias to decrease at larger 0^,(2.) (i.e. bias between AATS and CARL reduces to 10% 
if only CTep( 355 nm) > 0.15 Km' 1 are considered [Ferrare et al., this issue]). An unnoticed loss of 
sensitivity of the Raman lidar had occurred leading up to AIOP, and we expect bettor agreement 
from the recently restored system. However, the present comparison between AATS- 14 and 
CARL is valuable as it assesses the daytime retrievals of 0^(355 nm) of a Raman lidar in an 
operational setting. CARL is the only Raman lidar in the world designed to autonomously 
provide a continuous day and nighttime 10-year data record [Turner et al., 2002]. 

We emphasize the assessment of the uncertainties in the AATS-14 retrieved quantities. The 
instrument was calibrated at the Mauna Loa Observatory (MLO), Hawaii, 1.5 months before and 
1.5 months after the AIOP campaign and the calibration constants were fine tuned within the 
bounds of pre- and post-mission calibration by inspecting low t p (A.) spectra obtained near 
maximum flight altitude (5.6 km). 

The in-flight performance of AATS-14 was assessed by comparing r p (/ 1) obtained during 

low level legs (~90 m) against ground-based Sun photometers. The level of agreement between 
AATS-14 and two AERONET Sun photometers is similar to what we found from low altitude 
fly-bys over AERONET sites in previous campaigns. The agreement with a third ground-based 
Sun photometer (N1MFR) is particularly good, in fact operating four Sun photometers (including 
AATS-6) side-by-side on the ground in previous ARM IOPs did not result in a higher level of 
agreement. 

The robustness of the AATS-14 differentiated profiles of p w and aep(k) were tested by 
comparing p w obtained from AATS-14 and the aircraft in situ sensors. This presumes that an 
aircraft in-situ measurement of p w is more straightforward than measuring ambient a ep(X). 
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Averaged over 35 vertical profiles we find a relative rms difference of 20% and a small bias 
difference of 5% (in situ p* - AATS-14 p w ). 

Because most of the errors in AATS-14 t p (A.) or CWV are of systematic nature, they cancel 
out when differences (such as layer t p (k) or LWV) or differentiations (cr ep (A.) or p w ) are used. 
However, uncertainties in the AATS-14 cr ep (X) profiles arise from horizontal and temporal 
variability in the overlying aerosol. For the 25 o rp (X) vertical profiles used here this resulted in 
average Sa^JL = 453, 519, 675, 1558 nm)= 0.032, 0.029, 0.024, 0.014 Km' 1 equivalent to 101- 

176% if expressed as 8a ep (k)/a ep (k): The rms differences (other methods vs. AATS-14) in Table 
3 are smaller than these uncertainties. One might therefore conclude that the CepQC) 
measurements agree within the total error bars of AATS-14 alone. However, it is important to 
note, that the AATS-14 5a ep (^.) represent random errors. It is possible to turn these random 
errors into bias errors for a single profile by flying a ramped ascent or descent under a 
pronounced aerosol gradient without changing heading. However, none of the profiles used in 
this study was acquired using such a flight pattern. Moreover, averaged over an ensemble of 
profiles (as done in this study) we can rule out a systematic bias due to spatial variability as this 
would require flying each profile with the pattern described above, under gradients with the same 
mathematical sign. The same discussion applies to the uncertainties of AATS-14 layer x p (A.) 
where, in fact, we observe differences that are often larger than the random errors (see Figure 
10 ). 

Hence, we believe the observed biases to be statistically significant Furthermore, we find 
similar biases in the results published from previous field campaigns since 1996 involving 
AATS-6 or AATS-14: 
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Combining the results from AIOP with those from 5 previous field campaigns, we find 
airborne Nephelometer+PSAP measurements of layer r p (A = 450 - 700nm) to be biased slightly 

low (5 - 33%, average of 17%) when compared to airborne Sun photometer (AATS-6 or -14) 
values. 

From three previous campaigns we find layer x p (A.) calculated from airborne measurements 
of particle size distributions to be less than the AATS-6 or -14 values (average of 18%). 
However the data set for this computationally involved comparison is relatively small. 

In 5 previous field campaigns we have compared Oep(k) and t p (A.) vertical profiles from one 
of the AATS instruments with surface based or airborne lidars. Many comparisons result in small 
or no biases, however the biases that do occur are positive (i.e. lidar <jq,(X) larger than AATS 
values). 

There is a clear tendency for the remote sensing methods, lidar and airborne Sun 
photometers, to yield larger a ep (k) and i p (k) values than the in-situ methods. Adding the not 
previously used airborne cavity-ring-down technique (Cadenza instrument), did not significandy 
alter that tendency. In fact, Cadenza’s ability to in-situ measure at X=1550 nm, highlights the 
spectral signature of the low bias (i.e. low bias with respect to AATS- 14 increases with 
increasing wavelength). The low bias could be caused by particle sampling losses or incomplete 
corrections for shrinkage by evaporation of water, organics, or nitrates. 

We cite numerous studies, with no AATS involved, that also found the remote sensing 
methods (lidar and Sun photometers) to yield larger <7^(7.) or t p (k) values than the in-situ 
methods. 

Unknown gaseous absorption in the atmosphere (as postulated by Halthore et al. [1998]), 
not accounted for in the analysis of the Sun photometer and lidar data, could also lead to the 
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observed biases. However, Mlawer et al. [2000] provide “strong evidence that in this spectral 
range [350-1000 nm] there are no unmodeled molecular absorbers of significance to the 
atmospheric energy balance.” 

While we find that each of the methods investigated here has its strengths and weaknesses, 
there is no definitive proof that one of the methods is fun damentall y flawed. From the biases 
found in AIOP and previous studies, we conclude that the systematic error associated with 
measuring the tropospheric vertical profile of the ambient aerosol extinction with current state- 
of-the art instrumentation is 15-20% at visible wavelengths and potentially larger in the UV and 
near-infrared. Random errors, as measured by rms differences (e.g. Table 3), are considerably 
larger, ranging from 26% to 98%. 

Acknowledgments: The Atmospheric Radiation Measurement (ARM) Prog ram is 
sponsored by the U.S. Department of Energy, Office of Science, Office of Biological and 
Environmental Research, Environmental Sciences Division. The success of the Aerosol IOP was 
due to the hard work and dedicated efforts from a large team of scientists and investigators from 
National Laboratories and Universities, CIRPAS Twin Otter and Cessna pilots, crew and support 
personnel, SGP site personnel, ARM infrastructure support, weather forecaster, and support from 
Greenwood aviation at Ponca City airport. We thank ARM for the support of this IOP. This 
research was also supported by a grant form the NOAA Office of Global Programs. EJW and 
MPLNET are supported by the NASA Radiation Sciences Program and Earth Observing System. 
The authors would like to thank Jack Ji for operating AERONET and MPLNET instrumentation. 


33 


6 References 


Ackerman, T. and G. Stokes. 2003. "The Atmospheric Radiation Measurement Program." 
Physics Today 56: 38 - 45. 

Ackerman, A. S., O. B. Toon, D. E. Stevens, A. J. Heymsfield, V. Ramanathan, E. J. Welton, 
Reduction of tropical cloudiness by soot, Science, 288, 1042-1047, 2000. 

Ackermann, J., The extinction-to-backscatter ratio of tropospheric aerosol: A numerical study, J. 
Atmos. Oceanic Technol., 15, 1043-1050, 1998. 

Anderson, T. L. and J. A. Ogren, Determining Aerosol Radiative Properties Using the TSI 3563 
Integrating Nephelometer. Aerosol Science and Technology, 29, 57-69 (1998). 

Andrews E., P. J. Sheridan, J. A. Ogren, and R. Feirare, In Situ Aerosol Profiles over the 
Southern Great Plains CART site. Part I: Aerosol Optical Properties, J. Geophys. Res., 109, 
D06208, doi: 10. 1029/2003JD004025, 2004. 

Andrews E., P.J. Sheridan, M. Fiebig, A. McComiskey, J. A. Ogren, MJ. Bartholomew, D. 
Covert, R. Elleman, R. Gaspanni, D. Collins, H. Jonsson, B. Schmid, J. Wang, Comparison 
of methods for deriving aerosol asymmetry parameter. J. Geophys. Res. 2005 (this issue). 

Ansmann A., M. Riebesell, and C. Weitkamp. Measurement of atmospheric aerosol extinction 
profiles with a Raman lidar. OPTICS LETTERS, Vol. 15, No. 13 / July 1, 1990. 

Ansmann, A., et al.. Long-range transport of Saharan dust to northern Europe: The 11-16 
October 2001 outbreak observed with EARLINET, J. Geophys. Res., 108(D24), 4783, 
doi: 1 0. 1029/2003JD003757, 2003. 

Bane, J. M., Bluth, R., Flagg, C., Jonsson, H., Melville, W. K., Prince, M., Riemer, D. UNOLS 
now Oversees Research Aircraft Facilities for Ocean Science. Eos Trans. AGU, 85(41), 402, 
2004. 


34 



Bluth, R. T, Durkee, P. A., Seinfeld, J. H., Flagan, R. C., Russell, L. M., Crowley, P. A. and 
Finn, P. Center for Interdisciplinary Remotely-Piloted Aircraft Studies (CIRPAS). Bull. Amer. 
Meteor. Soc. 77, 2691-2699, 1996. 

Bogel, W. Neue Naherungsgleichungen fur den Sattigungsdruck des Wasserdampfes und fur die 
in der Meteorologie gebrauchlichen Luftfeuchte-Parameter. DLR-FB 77-52, Deutsche 
Forschungs imd Versuchsanstalt furLuft- und Raumfahrt, Oberpfaffenhofen (1977). 

Bond T. C., T. L. Anderson, and D. Campbell, Calibration and Intercomparison of Filter-Based 
Measurements of Visible Light Absorption by Aerosols. Aerosol Science and Technology, 
30:582-600(1999). 

Bucholtz A. et al.. Properties and Effects of Asian Aerosols over the Central California Coast 
During the ADAM-2003 (Asian Dust Above Monterey-2003) Field Study. A AAR 2003 
Annual Conference, Anaheim, California, October 20-24, 2003. 

Busch K.W. and M.A. Busch, (1999). “Cavity Ringdown Spectroscopy”, American Chemical 
Society, Washington, D.C. 

Campbell, J.R., D.L. Hlavka, E.J. Welton, C.J. Flynn, D.D. Turner, J.D. Spinhime, V.S. Scott, 
and IJL Hwang, 2002: Full-time, eye-safe cloud and aerosol lidar observations at 
Atmospheric Radiation Measurement program sites: instruments and data processing. J. 
Atmos. Oceanic Technol., 19, 431-442. 

Chameides, W. L., and M. Bergin, Soot takes center stage. Science, 297, 2214-2215, 2002. 

Chu, D. A., Y. J. Kaufman, G. Zibordi, J. D. Chem, J. Mao, C. Li, and B. N. Holben, Global 
monitoring of air pollution over land from the Earth Observing System-Terra Moderate 
Resolution Imaging Spectroradiometer (MODIS), J. Geophys. Res., 108(D21), 4661, 
doi:10. 1029/2002JD003 179, 2003. 


35 



Clough S. A., and M. J. Iacono, “Line-by-line calculations of atmospheric fluxes and cooling 
rates II: Application to carbon dioxide, ozone, methane, nitrous oxide, and the halocarbons,” 
J. Geophys. Res., 100, 16,519-16,535, 1995. 

Colarco P. et al. "Elevated injection height, long-range transport, and evolution of a Siberian 
forest fire smoke plume" to be submitted to J. Geophys. Res. (this issue). 

Collins, D. R., H. H. Jonsson, J. H. Seinfeld, R.C. Flagan, S. Gasso, D. A. Hegg, B. Schmid, P. 
B. Russell, J. M. Livingston, E. Ostrom, K. J. Noone, L. M. Russell, and J. P. Putaud, In situ 
aerosol size distributions and clear column radiative closure during ACE-2. Tellus, B 52, 498- 
525, 2000. 

Damoah, R., N. Spichtinger, C. Forster, P. James, I. Mattis, U. Wandinger, S. Beirle, T. Wagner, 
and A. Stohl, Around the world in 17 days — hemispheric-scale transport of forest fire smoke 
from Russia in May 2003. Atmospheric Chemistry and Physics, Vol. 4, pp 1311-1321, 23-8- 
2004. 

Delene, D.J. and J.A. Ogren. Variability of aerosol optical properties at four North American 
surface monitoring sites. J. Atmos. Sci., 59: 1 135-1150, 2002. 

Diner, D.J. WA. Abdou, J.E. Conel, K.A. Crean, B.J. Gaitley, M. Helmlinger, R.A. Kahn, J.V. 
Martonchik, S.H. Pilorz, and B.N. Holben, MISR aerosol retrievals over southern Africa 
during the SAFARI-2000 dry season campaign, Geophys. Res. Lett. 28, 3127-3130, 2001. 

Eck, T. F., B. N. Holben, J. S. Reid, O. Dubovik, A. Smirnov, N. T. O'Neill, I. Slutsker, and S. 
Kinne, Wavelength dependence of the optical depth of biomass burning, urban, and desert 
dust aerosols, J. Geophys. Res., 104, 31,333-31,349, 1999. 

Ferrare, R., S. Ismail, E. Browell, V. Brackett, M. Clayton, S. Kooi, S. H. Melfi, D. Whiteman, 
G. Schwemmer, K. Evans, P. Russell, J. Livingston, B. Schmid, B. Holben, L. Remer, A. 


36 


Smirnov, P. Hobbs. Comparisons of aerosol optical properties and water vapor among ground 
and airborne lidars and sun photometers during TARFOX. J. Geophys. Res., 105(D8), 9917- 
9933, 2000. 

Ferrare, R. A., D. D. Turner, L. A. Heilman, W. Feltz, O. Dubovik, and T. Tooman: Raman 
Lidar Measurements of the Aerosol Extinction-to-Backscatter Ratio Over the Southern Great 
Plains, J. Geophys. Res., 106, 20333-20347, 2001. 

Ferrare, R.A., Turner, D.D., Clayton, M„ Brackett, V., Tooman, T.T., Goldsmith, J.E.M., Ogren, 
J.A., Andrews, E„ Welton, EJ., Campbell, J.R., and Chin, M„ Vertical Variability of 
Aerosols and Water Vapor over the Southern Great Plains, DOE ARM Science Team 
Meeting, Broomfield, CO, 2003 

http://www.ar m. gov/ docs/documents/technical/conf 0304/fftrrare-ra.ndf 
Ferrare, R., D. Turner, M. Clayton, B. Schmid, J. Redemann, D. Covert, R. Elleman, J. Ogren, E. 
Andrews, J. Goldsmith, H. Jonsson, Raman Lidar Measurements of Aerosols and Water 
Vapor over the Southern Great Plains during the May 2003 Aerosol IOP. J. Geophys. 
Res., 2005 (this issue). 

Gao S„ D. A. Hegg, H. Jonsson, Aerosol Chemistry, and Light-Scattering and Hygroscopicity 
Budgets during Outflow from East Asia, Journal of Atmospheric Chemistry, Volume 46, 
Issue 1, September 2003, Pages 55 - 88. 

Gasso, S., D. A. Hegg, D. S. Covert, D. Collins, K. J. Noone, E. Ostrom, B. Schmid, P. B. 
Russell, J. M. Livingston, P. A. Durkee, and H. Jonsson, Influence of humidity on the aerosol 
scattering coefficient and its effect on the upwelling radiance during ACE-2. Tellus, B 52, 
546-567, 2000. 


37 



Hallar et al. ARM Aerosol IOP / IAP Flight Leg Comparisons -A Closure Experiment of 
Aerosol Optical Properties", to be submitted to JGR (this issue) 

Halthore R.N., S. Nemesure, S.E. Schwartz, D.G. Imre, A. Berk, E.G. Dutton, and M.H. Bergin, 
Models overestimate diffuse clear-sky surface irradiance: A case for excess atmospheric 
absorption”, Geophys. Res. Lett., Vol. 25, No. 19, 3591-3594, October 1, 1998. 

Hansen J., M. Sato, R. Ruedy, A. Lacis, and V. Oinas, Global warming in the twenty-first 
century: An alternative scenario. Proc. Natl. Acad. Sci., 97, 9875-9880, 2000. 

Harder J. W., J. W. Brault, P. V. Johnston, and G. H. Mount, Temperature dependent NO 2 cross 
sections at high spectral resolution, J. Geophys. Res., 102, 3861-3879, 1997. 

Hamson L., J. Michalsky, and J. Bemdt, 1994: Automated multifilter rotating shadow -band 
radiometer an instrument for optical depth and radiation measurements. Applied Optics, Vol. 
33, No. 22, 5118-5125. 

Hartley, W. S., Hobbs, P. V., Ross, J. L„ Russell, P. B., and Livingston, J. M. Properties of 
aerosols aloft relevant to direct radiative forcing off the mid-Atlantic coast of the United 
States. J. Geophys. Res., Vol. 105 , No. D8 , 9859-9886 (2000). 

Hegg D.A., J. Livingston, P .V. Hobbs, T. Novakov, and P. Russell, “Chemical apportionment of 
aerosol column optical depth off the mid-Atlantic coast of the United States”, J. Geophys. 
Res., Vol. 102, No. D21, 25’293-25’303, November 20, 1997. 

Holben, B. N., T. Eck, I. Slutsker, D. Tanre, J. B. Buis, A. Setzer, E. Vennote, J. A. Reagan, Y. 

J. Kaufman, T. Nakajima, F. Lavenu, I. Jankowiak, and A.Smimov, 1998: AERONET - A 
Federated Instrument Network and Data Archive for Aerosol Characterization. Rem. Sens. 
Env., 66, pp. 1-16. 


38 




Holben B. N. et al., An emerging ground-based aerosol climatology: Aerosol Optical Depth from 
AERONET, J. Geophys. Res., 106, 12067-12097, 2001. 

Jacobson, M. Z., Control of fossil-fiiel particulate black carbon and organic matter, possibly the 
most effective method of slowing global wanning, J. Geophys. Res., 107(D19), 4410, 
doi: 10. 1029/2001 JD001376, 2002. 

Kasten, F., Visibility in the phase of pre-condensation. Tellus, 21, 631-635, 1969. 

Kato, S., Bergin, M. H., Laulainen, N., Ferrare, R., Turner, D., Michalsky, Charlock, T. P., 
Clothiaux, E. E., Mace, G. G. and Ackerman, T. P. A comparison of the aerosol optical 
thickness derived from ground-based and airborne measurements. J. Geophys. Res., Vol 105, 
Dll 14,701-14,717,2000. 

Kaufman, Y., D. Tanre, and O. Boucher, A satellite view of aerosols in the climate system. 
Nature, 419, 215-223, 2002. 

Kaufman Y. J., J. M. Haywood, P. V. Hobbs, W. Hart, R. Kleidman and B. Schmid. Remote 
sensing of vertical distributions of smoke aerosol off the coast of Africa during SAFARI 
2000. Geophys. Res. Lett., Vol. 30, No. 16, doi:10.1029/2003GL017068, 2003. 

King M. D. and D. M. Byrne, 1976: A Method for Inferring Total Ozone Content from the 
Spectral Variation of Total Optical Depth Obtained with a Solar Radiometer. J. Atmos. Sci., 
Vol. 33, 2242-2251. 

Kotehenruther R.A., P.V. Hobbs, and D.A. Hegg, “Humidification factors for atmospheric 
aerosols off the mid-atlantic coast of the United States”, J. Geophys. Res. Vol 104, D2, 2239- 
2251,1999. 

Lelieveld, J., et al.. Global air pollution crossroads over the Mediterranean, Science, 298, 794- 
799, 2002. 


39 


Livingston, J. M., V. Kapustin, B. Schmid, P. B. Russell, P. K. Quinn, T. S. Bates, P. A. Durkee, 
P. J. Smith, V. Freudenthaler, D. S. Covert, S. Gasso, D. A. Hegg, D. R. Collins, R. C. 
Flagan, J. H. Seinfeld, V. Vitale, and C. Tomasi, Shipboard Sun photometer measurements of 
aerosol optical depth spectra and columnar water vapor during ACE 2 and comparison to 
selected land, ship, aircraft, and satellite measurements. Tellus, B 52, 594-619, 2000. 
Livingston J. M„ P. B. Russell, J. S. Reid, J. Redemann, B. Schmid, D. A. Allen, O. Torres, R. 
C. Levy, L. A. Remer, B. N. Holben, A. Smirnov, O. Dubovik, E. J. Welton, J. R. Campbell, 
J. Wang, S. A. Christopher, Airborne Sun photometer measurements of aerosol optical depth 
and columnar water vapor during the Puerto Rico Dust Experiment, and comparison with 
land, aircraft, and satellite measurements, J. Geophys. Res., 108 (D19), 8588, 
doi : 1 0. 1029/2002JD002520, 2003. 

Livingston J. M., B. Schmid, P. B. Russel], J. A. Eilers, R. W. Kolyer, J. Redemann, S. A. 
Ramirez, J-H. Yee, W. H. Swartz, C. R. Trepte, L. W. Thomason, M. C. Pitts, M. A. Avery, 
C. E. Randall, J. D. Lumpe, R. M. Bevilacqua, M. Bittner, T. Erbertseder, R. D. McPeters, R. 
E. Shetter, E. V. Browell, J. B. Kerr, K. Lamb. Retrieval of ozone col umn content from 
airborne Sun photometer measurements during SOLVE II: Comparison with coincident 
satellite and aircraft measurements. Atmospheric Chemistry and Physics Discussions, 5, 243- 
286, 2005. 

Magi, B. L, P. V. Hobbs, B. Schmid, and J. Redemann, Vertical profiles of light scattering, light 
absorption and single scattering albedo during the dry, biomass burning season in southern 
Africa and comparisons of in situ and remote sensing measurements of aerosol optical depths, 

J. Geophys. Res., 108(D13), 8504, doi:10.1029/2002JD002361, 2003. 


40 


Magi B. I., P. V. Hobbs, T. W. Kirchstetter, T. Novakov, D. A. Hegg, S. Gao, J. Redexnann, B. 
Schmid, Aerosol Properties and Chemical Apportionment of Aerosol Optical Depth at 
Locations off the United States East Coast in July and August 2001. J. Atmos. Sci., in press, 
2005. 

Masonis S J, K. Franke, A. Ansmann, D. Muller, D. Althausen, J.A. Ogren, A. Jefferson, P. J. 
Sheridan, An intercomparison of aerosol light extinction and 180°-backscatter as derived 
using in situ instruments and Raman lidar during the INDOEX field campaign, J. Geophys. 
Res., 107, No. D19, 8014, doi:10.1029/2000JD000035, 2002. 

Matsumoto, T., P. B. Russell, C. Mina, W. Van Ark and V. Banta, 1987: Airborne Tracking Sun 
photometer. J. Atmos. Ocean. Tech., Vol. 4, 336-339. 

McGill M., D. Hlavka, W. Hart, J. Spinhime, S. Scott, and B. Schmid. The Cloud Physics Lidar: 
Instrument Description and Initial Measurement Results. Applied Optics- LP, Volume 41, 
Issue 18, 3725-3734, June 2002. 

Menon, S., J. Hansen, L. Nazarenko and Y. Luo, Climate effects of black carbon aerosols in 
China and India, Science, 297, 2250-2252, 2002. 

Mlawer E. J., P. D. Brown, S. A. Clough, L. C. Harrison, J. J. Michalsky, P. W. Kiedron, and T. 
Shipper!, Comparison of spectral direct and diffuse solar irradiance measurements and 
calculations for cloud-free conditions, Geophys. Res. Lett.. 27, 2,653-2,656, 2000. 

Murayama T., S.J. Masonis, J. Redemann, T.L. Anderson, B. Schmid, J.M. Livingston, P.B. 
Russell, B. Huebert, S.G. Howell, C.S. McNaughton, A. Clarke, M. Abo, A. Shimizu, N. 
Sugimoto, M. Yabuki, H. Kuze, S. Fukagawa, K.L. Maxwell, R.J. Weber, D.A. Orsini, B. 
Blomquist, A. Bandy, D. Thornton, An intercomparison of lidar-derived aerosol optical 


properties with airborne measurements near Tokyo during ACE-Asia. J. Geophys. Res., 
108(D23), 8651, doi:10.1029/2002JD003259, 2003. 

O'Keefe A. and Deacon D.A.G. Cavity ring-down optical spectrometer for absorption 
measurements using pulsed laser sources. Rev. Sci. Inst., vol. 59, pg- 2544 (1988). 

Penner, J. E., C omment on “Control of fossil-fuel particulate black carbon and organic matter, 
possibly the most effective method of slowing global warming” by M. Z. Jacobson, J. 
Geophys. Res., 108(D24), 4771, doi:10.1029/2002JD003364, 2003. 

Penner, J. E., S. Y. Zhang, and C. C. Chuang. Soot and smoke aerosol may not warm climate. J. 

Geophys. /to.,108(D21),4657, 10. 1029/2003 JD003409, 2003. 

Petzold, A., M. Fiebig, H. Flentje, A. Keil, U. Leiterer, F. Schroder, A. Stifter, M. Wendisch, and 
P. Wendling (2002), Vertical variability of aerosol properties observed at a continental site 
during the Lindenberg Aerosol Characterization Experiment (LACE 98), J. Geophys. Res., 
107(021), 8128, doi:10.1029/2001JD001043. 

Ramanathan, V., PJ. Crutzen, J.T. Kiehl, and D. Rosenfeld, Aerosol, Climate and the 

Hydrological Cycle. Science, 294, 2119-2124, 2001. 

Redemann, J., R. P. Turco, K. N. Liou, P. B. Russell, R. W. Bergstrom, B. Schmid, J. M. 
Livingston, P. V. Hobbs, W. S. Hartley, S. Ismail, R. A Ferrare, E. V. Browell, Retrieving the 
Vertical Structure of the Effective Aerosol Complex Index of Refraction From a Combination 
of Aerosol In Situ and Remote Sensing Measurements During TARFOX. J. Geophys. Res., 
105(08), 9949-9970, 2000 

Redemann J., P.B. Russell, and P. Hamill, Dependence of aerosol light absorption and single- 
scattering albedo on ambient relative humidity for sulfate aerosols with black carbon cores, J. 
Geophys. Res., 106, 27,485-27,495, 2001. 


42 


Redemann, J., S. J. Masonis, B. Schmid, T. L. Anderson, P. B. Russell, J. M. Livingston, O. 
Dubovik, and A. D. Clarke, Cl ear-column closure studies of aerosols and water vapor aboard 
the NCAR C-130 during ACE-Asia, 2001, J. Geophys. Res., 108{ D23), 8655, 
doi: 10. 1029/2003JD003442, 2003. 

Redemann, J., B. Schmid, J.A. Eilers, R. Kahn, R.C. Levy, P.B. Russell, J.M. Livingston, P.V. 
Hobbs, W.L. Smith Jr., B.N. Holben, Suborbital measurements of spectral aerosol optical 
depth and its variability at sub-satellite-grid scales in support of CLAMS, 2001, J. Atmos. 
Sci., in press 2005. 

Remer L.A., S. Gasso, D.A.Hegg, YJ. Kaufman, and B.N. Holben, “Urban/industrial aerosol: 
Ground-based Sun/sky radiometer and airborne in situ measurements”, J. Geophys. Res., Vol. 
102, No. D14, 16’ 849-16’ 859, July 27, 1997. 

Rothman L.S., K. Chance, J. Schroeder, and A. Goldman. New Edition of HTTRAN Database. 

1 1 th ARM Science Team Meeting Proceedings, Atlanta, Georgia, March 19-23, 2001. 
Rothman L.S., and J. Schroeder, Millenium HlTKAN Compilation. 12 th ARM Science Team 
Meeting Proceedings, St. Petersburg, Florida, April 8-12, 2002. 

Rotstayn, L. D., and U. Lohmann, Tropical rainfall trends and the indirect aerosol effect, J. 
Clim., 15, 2103-2116, 2002. 

Russell, P. B., J. M. Livingston, E. G. Dutton, R. F. Pueschel, J. A. Reagan, T. E. Defoor, M. A. 
Box, D. Allen, P. Pilewskie, B. M. Herman, S. A. Kinne, and D. J. Hofmann, Pinatubo and 
pre-Pinatubo optical-depth spectra: Mauna Loa measurements, comparisons, inferred particle 
size distributions, radiative effects, and relationship to lidar data. J. Geophys. Res., 98, 
22'969-22'985, 1993a. 


43 


Russell, P.B., J. M. Livingston, R. F. Pueschel, J. A. Reagan, E.V. Browell, G. C. Toon, P.A. 
Newman, M.R. Schoeberl, L.R. Lait, L. Pfister, Q. Gao, and B. M. Herman, "Post-Pinatubo 
Optical Depth Spectra vs. Latitude and Vortex Structure: Airborne Tracking Sun photometer 
Measurements in AASE II," Geophys. Res. Lett., 20, 2571-2574, 1993b. 

Russell, P. B., P. V. Hobbs, and L. L. Stowe, Aerosol properties and radiative effects in the 
United States Mid-Atlantic haze plume: An overview of the Tropospheric Aerosol Radiative 
Forcing Observational Experiment (TARFOX), J. Geophys. Res., 104, 2213-2222, 1999a. 

Russell, P. B., J. M. Livingston, P. Hignett, S. Kinne, J. Wong, and P. V. Hobbs, Aerosol- 
induced radiative flux changes off the United States Mid- Atlantic coast. Comparison of values 
calculated from Sun photometer and in situ data with those measured by airborne 
pyranometer, J. Geophys. Res., 104, 2289-2307, 1999b. 

Russell, P. B., J. M. Livingston, O. Dubovik, S. A. Ramirez, J. Wang, J. Redemann, B. Schmid, 
M. Box, and B. N. Holben, Sunlight transmission through desert dust and marine aerosols: 
Diffuse light corrections to Sun photometry and pyrheliometry, J. Geophys. Res., 109, 
D08207, doi: 10. 1 029/2003 JD004292, 2004a 

Russell P., J. Livingston, B. Schmid, J. Eilers, R. Kolyer, J. Redemann, S. Ramirez, J-H. Yee, W. 
Swartz, R. Shetter, C. Trepte, A. Risley, Jr., B. Wenny, J. Zawodny, W. Chu, M. Pitts, J. 
Lumpe, M. Fro mm ,C. Randall, K. Hoppel, R. Bevilacqua. "Aerosol optical depth 
measurements by airborne Sun photometer in SOLVE II: Comparisons to SAGE HI, POAM 
HI and airborne spectrometer measurements" Atmospheric Chemistry and Physics 
Discussions, 4, 7291-7353, 2004b 


44 



Sato M., J. Hansen, D. Koch, A. Lacis, R. Ruedy, O. Dubovik, B. Holben, M. Chin, and T. 
Novakov, Global atmospheric black carbon inferred from AERONET. Proc. Natl. Acad. Sci., 
100, 6319-6324, 2003. 

Schmid, B., and C. Wehrli, Comparison of sun photometer calibration by Langley technique and 
standard lamp, Appl. Opt., 34, 4500-4512, 1995. 

Schmid, B., K. J. Thome, P. Demoulin, R. Peter, C. Matzler, and J. Sekler, Comparison of 
modeled and empirical approaches for retrieving columnar water vapor from solar 
transmittance measurements in the 0.94 micron region, J. Geophys. Res., 101, 9345-9358, 
1996. 

Schmid, B., C. Matzler, A. Heimo, and N. Kampfer, Retrieval of Optical Depth and Size 
Distribution of Tropospheric and Stratospheric Aerosols by Means of Sun Photometry. IEEE 
Geosci. Remote. Sens., 35(1), 172-182, 1997. 

Schmid, B., P. R. Spyak, S. F. Biggar, C. Wehrli, J. Sekler, T. Ingold, C. Matzler, and N. 
Kampfer, Evaluation of the applicability of solar and lamp radiometric calibrations of a 
precision Sun photometer operating between 300 and 1025 nm, Appl. Opt., 37, 3923-3941, 
1998. 

Schmid B., J. Michalsky, R. Halthore, M. Beauhamois, L. Harrison, J. Livingston, P. Russell, B. 
Holben, T. Eck, and A. Smirnov, Comparison of Aerosol Optical Depth from Four Solar 
Radiometers During the Fall 1997 ARM Intensive Observation Period, Geophys. Res. Lett., 
26(17), 2725-2728, 1999. 

Schmid, B., J. M. Livingston, P. B. Russell, P. A. Durkee, H. H. Jonsson, D. R. Collins, R. C. 
Flagan, J. H. Seinfeld, S. Gasso, D. A. Hegg, E. Ostrom, K. J. Noone, E. J. Welton, K. J. 
Voss, H. R. Gordon, P. Formenti, and M. O. Andreae, Clear sky closure studies of lower 


45 



tropospheric aerosol and water vapor during ACE 2 using airborne Sun photometer, airborne 
in-situ, space-borne, and ground-based measurements, Tellus, B 52, 568-593, 2000. 

Schmid B., J. J. Michalsky, D. W. Slater, J. C. Barnard, R. N. Halthore, J. C. Liljegren, B. N. 
Holben, T. F. Eck, J. M. Livingston, P. B. Russell, T. Ingold, and I. Slutsker, Comparison of 
columnar water-vapor measurements from solar transmittance methods. Appl. Opt., 40, 1886- 
1896,2001. 

Schmid B., J. Redemann, P. B. Russell, P. V. Hobbs, D. L. Hlavka, M. J. McGill, B. N. Holben, 
E- J- Welton, J. R. Campbell, O. Torres, R. A. Kahn, D. J. Diner, M. C. Helmlinger, D. A. 
Chu, C. Robles Gonzalez, and G. de Leeuw, Coordinated airborne, spaceborne, and ground- 
based measurements of massive, thick aerosol layers during the dry season in southern Africa, 
J. Geophys. Res, 108(D13), 8496, doi:10.1029/2002JD002297, 2003a. 

Schmid B., D. A. Hegg, J. Wang, D. Bates, J. Redemann, P. B. Russell, J. M. Livingston, H. H. 
Jonsson, E. J. Welton, J. H. Seinfeld, R. C. Flagan, D. S. Covert, O. Dubovik, and A. 
Jefferson. Column closure studies of lower tropospheric aerosol and water vapor during ACE- 
Asia using airborne Sun photometer, airborne in-situ and ship-based lidar measurements. J. 
Geophys. Res., 70S(D23), 8656, doi:10.1029/2002JD003361, 2003b. 

Sheridan P. J., W. P. Amott, J. A. Ogren, E. Andrews, D. B. Atkinson, D. S. Covert, H. 
Moosmiiller, A. Petzold, B. Schmid, A. W. Strawa, R. Varma, A. Virkkula. The Reno Aerosol 
Optics Study: An Evaluation of Aerosol Absorption Measurement Methods. Aerosol Science 
& Technology, Volume 39, Number 1, January 2005, pp. 1-16, doi: 
10.1080/027868290901891. 

Sheridan P.J., D.J. Delene, and J.A. Ogren, Four years of continuous surface aerosol 
measurements from the Department of Energy’s Atmospheric radiation Measurement 


46 



Program Southern Great Plains Cloud and Radiation Testbed site. J. Geophys. Res., .106, D18, 
20735-20747, 2001. 

Sivaraman, C., DJD. Turner, and C.J. Flynn, 2004: Techniques and methods used to determine 
the aerosol best estimate value-added product at the SGP central facility. Proceedings of the 
14th ARM science team meeting, Albuquerque, NM. Available at 
http://www .arm. gov/publications/proceedings/conf 1 4/extended abs/sivaraman-c .pdf 

Smirnov, A., B. N. Holben, T. F. Eck, O. Dubovik, and I. Slutsker, Cloud-screening and quality 
control algorithms for the AERONET database, Rem. Sens. Environ., 73, 337- 349, 2000. 

Smirnov, A, AERONET Processing Algorithms Refinement, AERONET Workshop, May 10 - 
14, 2004, El Arenosillo, Spain. 

Smith J. D. and D. B. Atkinson, "A Portable Pulsed Cavity Ring-Down Transmissometer for 
Measurement of the Optical Extinction of the Atmospheric Aerosol”, Analyst, 126, 1216, 
( 2001 ). 

Spinhime, J.; Welton, E.; Palm, S.; Hlavka, D.; Hart, W.; Mahesh, A.; Atmospheric 
measurements by the geoscience laser altimeter system: initial results. Geoscience and 
Remote Sensing Symposium, 2003. IGARSS '03. Proceedings. 2003 IEEE International, 
Volume: 3, 21-25 July 2003, Pages: 1537 - 1539. 

Spinhime, J. D., J. A. R. Rail, and V. S. Scott, Compact Eye Safe Lidar Systems, Rev. Laser 
Eng., 23, 112-118, 1995. 

Sprent and Dolby (1980). The Geometric Mean Functional Relationship. Biometrics 36: 547- 
550. 


47 



Strawa, A. W., R. Castaneda, T. Owano, D. S. Baer, and B. A. Paldus, 2003: The measurement 
of aerosol optical properties using continuous wave cavity ring-down techniques. J. Atmos. 
Oceanic Tech., 20, 454-A65. 

Strawa, A.W., Amott, A.P., Elleman, R., Hallar, A.G., Covert, D., Ricci, K., Provencal, R., 
Owano,T., Jonsson, H, Schmid, B., Luu, A.P., Bokarius, K. In-Situ Measurement of Aerosol 
Extinction and Scattering Made During the DOE Aerosol IOP. J. Geophys. Res. 2005 (this 
issue). 

Torres, O., P. K. Bhartia, J. R. Herman, A. Sinyuk, P. Ginoux and B. Holben, A long-term record 
of aerosol optical depth from TOMS observations and comparison to AERONET 
measurements, J. Atmos. Sci., 59, 398-413, 2002. 

Turner, D.D., R.A. Ferrare, and L.A. Brasseur, Average aerosol extinction and water vapor 
profiles over the Southern Great Plains. Geophys. Res. Lett., 28, 4441-4444, 2001. 

Turner D.D., R.A. Ferrare, L.A. Heilman, W.F. Feltz, and T. Tooman, Automated Retrievals of 
Water Vapor and Aerosol Profiles over Oklahoma from an Operational Raman Lidar, J. 
Atmos. Oceanic Tech., 19, 37-50, 2002. 

VanCuren, R. A., Asian aerosols in North America: Extracting the chemical composition and 
mass concentration of the Asian continental aerosol plume from long-term aerosol records in 
the western United States, J. Geophys. Res., 108(D20), 4623, doi:10.1029/2003JD003459, 
2003. 

Virkkula A., N, C. Ahlquist, D. S. Covert, W. P. Amott, P. J. Sheridan, P. K. Quinn, D. J. 
Coffman, Modification, Calibration and a Field Test of an Instrument for Measuring Light 
Absorption by Particles. Aerosol Science & Technology, Volume 39, Number 1, January 
2005, pp. 68-83, doi: 10.1080/027868290901963. 


48 



Wang J., R. C. Flagan, J. H. Seinfeld, H. H. Jonsson, D. R. Collins, P. B. Russell, B. Schmid, J. 
Redemann, J. M. Livingston, S .Gao, D. A. Hegg, E. J. Welton, and D. Bates. Clear-column 
radiative closure during ACE-Asia: Comparison of multiwavelength extinction derived from 
particle size and composition with results from sunphotometry, J. Geophys. Res., 107(D23), 
4688, doi: 10. 1029/2002JD002465, 2002. 

Wang, J, Christopher, S A, Nair, U S, Reid, J, Prins, E M, Szykman, J. (2004). The Effect of 
Central American Smoke Aerosols on the Air Quality and Climate over the Southeastern 
United States: First Results from RAMS-AROMA. Eos Trans. AGU, 85(47), Fall Meet. 
Suppl., Abstract A1 1 A-0022. 

Welton, E.J., K.J. Voss, H.R. Gordon, H. Maring, A. Smirnov, B. Holben, B. Schmid, J.M. 
Livingston, P.B. Russell, P.A. Durkee, P. Formenti, M.O. Andreae, "Ground-based Lidar 
Measurements of Aerosols During ACE-2: Instrument Description, Results, and Comparisons 
with other Ground-based and Airborne Measurements", Tellus B, 52, 635-650, 2000. 

Welton, E. J., J. R. Campbell, J. D. Spinhime, and V. S. Scott, Global monitoring of clouds and 
aerosols using a network of micro-pulse lidar systems ", in Lidar Remote Sensing for Industry 
and Environmental Monitoring, U. N. Singh, T. Itabe, N. Sugimoto, (eds.), Proc. SPIE, 4153, 
151-158,2001. 

Welton, E.J., K.J. Voss, P.K. Quinn, P.J. Flatau, K. Markowicz, J.R. Campbell, J.D. Spinhime, 
H.R. Gordon, and J.E. Johnson, "Measurements of aerosol vertical profiles and optical 
properties during INDOEX 1999 using micro-pulse lidars", J. Geophys. Res., 107, 8019, 
doi: 10. 1029/2000 JD000038, 2002. 


49 


Table 1: Comparison of airborne and surface measured AODs during low altitude flybys. 
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Table 2: Comparison of water vapor retrieved from AATS-14 and EdgeTech 137-C3 chilled mirror. 














Table 4 : Comparison of layer x p (k) from six different methods 
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Table 6 Comparison of layer t p (A.) measured with AATS-14 (or AATS-6) and calculated from size distributions 
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7 Figure captions: 

Figure 1: Selection of t p (X) vertical profiles from AIOP. 

Figure 2: Vertical profiles of <r ep (X) derived from the x p (X) profiles shown in Figure 1. 

Figure 3: Vertical profiles of CWV for cases shown in Figure 1. 

Figure 4: Vertical profiles of p w from AATS-14 and EdgeTech 137-C3 chilled mirror sensor for cases 
shown in Figure 3. 

Figure 5: Comparison of p w from AATS-14 and EdgeTech 137-C3 chilled mirror sensor for 35 vertical 
profiles. 

Figure 6: Comparison of LWV from AATS-14 and EdgeTech 137-C3 chilled mirror sensor for 35 vertical 
profiles. AATS-14 error bars are based on horizontal distance spanned by a profile, combined with 
average horizontal variability of CWV in AIOP flights. 

Figure 7: Vertical profiles of a ep { 675 nm) from Nephelometer+PSAP and AATS-14 for the cases shown 
in Figure 2. (The Cadenza a ep (675nm) profiles are virtually indistinguishable from Nephelometer+PSAP 
data points and are therefore not plotted) 

Figure 8: Comparison of a ep (675 nm) from Nephelometer+PSAP and Cadenza for all 26 vertical profiles. 
Black: 1:1 line, blue: regular y vs. x regression, green: inverted x vs. y regression, red: bisector of blue 
and green lines (i.e. least squares bi-sector method [Sprent and Dolby, 1980]). 

Figure 9: Comparison of o ep (675 nm) from AATS-14 and Cadenza for all 26 vertical profiles. Regression 
lines as in Figure 8. 

Figure 10: Comparison of layer t p (675 nm) from AATS-14 and Cadenza for all 26 vertical profiles. 
Regression lines as in Figure 8. AATS-14 error bars are based on horizontal distance spanned by a profile, 
combined with average horizontal variability of AOD in AIOP flights. Cadenza error bars reflect 10% 
uncertainty. 
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Figure 11: Ratios of humidification factors (see equation (4)) from a 1-year analysis (March 2000 - 
February 2001) of surface-based dry / humidified a sp Q^ = 450, 550, 700 nm) (submicron particles only) 
measured with a TSI nephelometer at the SGP CRF. Solid line: X=700 nm / 550 nm, dashed line A.=450 
nm / 550 nm. 

Figure 12: Vertical profiles of a ep from Nephelometer+PSAP (519 nm), MPLNET (523 nm) and 
MPLARM (523 nm) for the cases shown in Figure 2. 

Figure 13: Vertical profiles of x p from MPLARM (523 nm), MPLNET (523 nm) and AATS-14 (519 nm) 
for the cases shown in Figure 2. 
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How weU can we measore the vertical profile of tropospheric aerosol extinction? 


Aiunors: 
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Michalsky , A, Smirnov, B. Holben, and J. Baraard 
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Apr ™? T^f ^ Department of Energy Atmospheric Radiation Measurement (ARM) 
Aerosol Intensive Operations Period (AIOP, May 2003) yielded one of the best 
measurement sets obtained to-date to assess our ability to measure the vertical profile of 
ambient aerosol exoncbon oJX, in the lower troposphere During one a heavily 

instrumented aircraft with well characterized aerosol sampling ability carrying well 
proven and new aerosol instrumentation, devoted most of the 60 available flight hours to 
flying vertical profiles over the heavily instrumented ARM Southern Great Phans (SGP) 

™ s “ “ 10 compare -SSSTrSffi 

obtained from 6 different instruments: airborne Sun photometer (AATS-14) airborne 

e,astlc backscatter lidars. We find the in-situ measured 
a v (X) to be lower than the AATS-14 denved values. Bias differences are 0 002 - 0 004 

™ 12 " 17% m thC Visible ’ 0r 45% in the near-inftaredL On the other hand, 

n nnl t ^ ^ t0 46 hdar °«A) are higher: Bias differences are 

’ ( 1 3 and 0.007 Km (24%) for the two elastic back-scatter lidars (MPLNET 

and MPLARR ^ =523 ^ and 0 029 Km (54%) for ^ lidar (X=355 nm) An 

^J° SS ° f SenSltlVlty of * e Raman lidar had occurred leading up to AIOP and we 
expect better agreement from the recently restored system. 

Looking at the collective results from 6 field campaigns conducted since 1996 
airborne m situ measurements of a#.) tend to be biased slightly low (17% at visible 
wavelengths ) when compared to airborne Sun photometer ojk). On the other hand, 
o„(A) values denved from lidars tend to have no or positive biases. 

we conc,ude that tbe typical systematic error associated 
P° SphenC vertical V° me ofthe ambient aerosol extinction with 

larger Ltoewtd^-tS^ “ 15 ' 2 ° % * “ d *"*** 
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